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Abstract. In this paper, both structural and dynamical stabilities of steady tran- 
sonic shock solutions for one-dimensional Eulcr-Poission system are investigated. First, 
a steady transonic shock solution with supersonic backgroumd charge is shown to be 
structurally stable with respect to small perturbations of the background charge, pro- 
vided that the electric field is positive at the shock location. Second, any steady transonic 
shock solution with the supersonic background charge is proved to be dynamically and 
exponentially stable with respect to small perturbation of the initial data, provided the 
electric field is not too negative at the shock location. The proof of the first stability 
result relies on a monotonicity argument for the shock position and the downstream 
density, and a stability analysis for subsonic and supersonic solutions. The dynamical 
stability of the steady transonic shock for the Euler-Poisson equations can be transformed 
to the global well-posedness of a free boundary problem for a quasilinear second order 
equation with nonlinear boundary conditions. The analysis for the associated linearized 
problem plays an essential role. 



1. Introduction and Main Results 
The following system of one dimensional Euler-Poisson equations: 

Pt + {pu)x = 0, 



(1) 



{pu)t + {p{p) + pu'^)x = pE, 
Ex = p- b{x), 



models several physical flows including the propagation of electrons in submicron semicon- 
ductor devices and plasmas (cf. |2H])( hydrodynamic model), and the biological transport 
of ions for channel proteins (cf. [2]). In the hydrodynamical model of semiconductor 
devices or plasma, u, p and p represent the macroscopic particle velocity, electron den- 
sity and pressure, respectively, E is the electric field, which is generated by the Coulomb 
force of particles. b{x) > stands for the density of fixed, positively charged background 

ions. The biological model describes the transport of ions between the extracellular side 
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and the cytoplasmic side of the membranes (|2]). In this case, p, pu and E are the ion 
concentration, the ions translational mass, and the electric field, respectively. 

In this paper, we prove two distinct stability results for steady transonic shocks These 
are solutions of the following time-independent equations 



The first result concerns the stability of such solutions under perturbation of a constant 
background charge density. That is a purely stationary result. The second concerns the 
global in time stability for solutions whose initial data are close to a stationary solution. 
We assume that p satisfies: 



First, consider the boundary value problem for in an interval < x < L with the 
boundary conditions: 



We assume ui > Q and Ur > 0. By the first equation in ([2]), we know that pu{x) = 
constant{0 < x < L), so the boundary data must satisfy 



(2) 



{pu)cc = 0, 

< (pip) + pu% = pE, 
E^ = p-h{x). 



(3) 



p(0) = p'(0) = 0, p'{p) > 0, p"{p) > 0, forp > 0, p{+oo) = +oo. 



(4) 



(p,M,E)(0) = {pi,Ul,Ei), {p,u){L) = {pr,Ur). 



(5) 



PlUl 



PrUr- 



If one denotes 



(6) 



PlUi = PrUr = J, 



then pu{x) = J(0 < x < L) and the velocity is given by 



(7) 



u = J/p. 



Thus the boundary value problem reduces to 




with the boundary conditions: 



(9) 



{p,Em = {pi,Ei), p{L)=pr. 



Use the terminology from gas dynamics to call c = \/p'{p) the sound speed. There is a 
unique solution p = Ps for the equation 

(10) p'ip) = .P/p\ 

which is the sonic state (recall that J = pu). Later on, the flow is called supersonic if 

(11) p'{p) < J^/p^, i.e. p<ps. 
Similarly, if 

(12) p\p) > .P/p\ I.e. p > Ps, 

then the flow is said to be subsonic. 

A piecewise smooth solution (p, E) with p > of (8) (or equivalently (2) with u = ^) 
is said to be a transonic shock solution, if it is separated by a shock discontinuity, and of 
the form 

{Psup,Esup){x), 0<X<Xo, 
{psub, Esub){x), Xo < X < L, 

satisfying the Rankine-Hugoniot conditions 

P{psup{Xo-)) H 7 T = P{psubiXo + )) H —,Esup{Xo-) = Esub{Xo + ), 

Psup[Xo — ) Psub[Xo + ) 

and is supersonic behind the shock and subsonic ahead of the shock, i.e., 

(13) PsupiXo-) < Ps < PsubiXo + ), 

and 

(14) psubix) > ps, for X e [xo,L]. 

(We assume that the velocity is always positive.) 

Note that (IHI)i is singular at the sonic state {p'{ps) — = 0) and the coefficient of 
Px changes sign for the supersonic flow and subsonic flow. This makes the problem of 
determining which kind of boundary conditions should be posed to make the boundary 
value problem well-posed a subtle one. In the previous works, some pure subsonic or 
supersonic solutions are obtained for both 1-dimensional and multidimensional cases (cf. 
[3 El [31]). For a viscous approximation of transonic solutions in 2-d case for the equations 
of semiconductors, see frU\ . However, there have been only a few results for the transonic 
flow. In the following, we list several results which are closely related to the present paper. 
First, a boundary value problem for ([8]) was discussed in [1] for a linear pressure function 
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of the form p{p) = kp, furthermore, the boundary conditions read p(0) = p{L) = p with 
p being a subsonic state and the density of the background charge satisfied < b < pg. 
The solution in jlj may contain transonic shock. On the other hand, since the boundary 
conditions and the pressure function are special in p.J, it is desired to consider the more 
general boundary conditions with more general equation of states. In |37], a phase plane 
analysis is given for system ([8]). However, no transonic shock solutions were constructed 
in [37]. A transonic solution which may contain transonic shocks was constructed by 
Gamba (cf. [9]) by using a vanishing viscosity limit method. However, the solutions as 
the limit of vanishing viscosity may contain boundary layers. Therefore, the question 
of well-posedness of the boundary value problem for the inviscid problem can not be 
answered by the vanishing viscosity method. Moreover, the structure of the solutions 
constructed by the vanishing viscosity method in [9J is shown to be of bounded total 
variation and possibly contain more than one transonic shock. 

A thorough study of the transonic shock solutions for one-dimensional Euler-Poisson 
equations with a constant doping profile b{x) = Bq was given in [25] for both cases when 
< bo < ps and bo > pg. In the present paper, we concentrate on the case when 
< 6o < Ps- In [25j, when < 6o < Ps, and if the boundary data {pi,Ei), pr and the 
interval length L satisfy some conditions (see p5] for details), then the problem ([8]) and 
([9]) admits a unique transonic shock solution. The stability of the transonic solutions 
obtained in [23] when 6 is a small perturbation of given background charge bo = bo{x), is 
analyzed in the following theorem. 

Theorem 1. Let J > be a constant, and let bo satisfy 
(15) < min bo(x) < max bo(x) < p. 

and {pi,Ei) be a supersonic state (0 < pi < ps), Pr be a subsonic state (pr > Ps)- If the 
boundary value problem ^ and ^ admits a unique transonic shock solution {p^'^\ E^^^) 
for the case when b{x) = bo{x) (x G [0,L]j with a single transonic shock located at x = 
xo G (0, L) satisfying 



(16) 



EW(xo+) = ^(°Ha;o-)>0, 



then there exists cq > such that if 



(17) 



b - &o||co[o,L] = e < eo, 
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then the boundary problem ^ and ^ admits a unique transonic shock solution (p, E) 
with a single transonic shock locating at some Xq G [xq — Ce,XQ + Ce] for some constant 
C>0. 

Remark 1. When Bq^x) = const, it should be noted that there are a large class of boundary 
data which ensure the existence and uniqueness of the transonic shock solutions satisfying 
the assumptions in TheoremUl , see [25] . 

Second, we would like to investigate the dynamical stability of the steady transonic 
shock solutions. For a given function b{x) satisfying < 6 = b{x) < ps for x G [0, L], and 
a constant J > 0, let 

{p-,J/p-,E^){x), if < X < xo, 
(p+, J/p+,E+)(x), if xo < X < L 

be a steady transonic shock solution of ([2]) satisfying the boundary conditions 
(19) {p,Em = {pi,Ei), p{L) = p,. 

Precisely, we assume that {p,u,E) is supersonic for < x < Xq, and subsonic for xq < 
X < L, i.e., 



{p,u,E){x) 



(20) 



p'{p-){x) < J^/p^ (x), for < X < Xq, 
p'{p+){x) > P/p^(x), for Xq < X < L. 

At X = Xo, the Rankine-Hogoniot conditions are satisfied: 



(21) ( p(p_) + —] (xo) = (pip+) + — ] (xo),^-(xo) = E+ixo). 



P-J \ Ph 

Finally, we assume the solution is away from vacuum: 

(22) inf p(x) > 0. 

xg[0,L] 

It is easy to see that we can extend {p_,E_) to be a smooth supersonic solution of ([2]) 
on [0,Xo + S] for some 6 > 0, which coincides with {p^^EJ) on [0,Xo]. Later on, we still 
use {p^.EJ) to denote this extended solution. Similarly, we will denote {p+,E^) to be a 
subsonic solution of ([2]) on [xq — S,L] for some 6 > 0, which coincides with {p+,E+) in 
dUD on [xo,L]. 

We consider the initial boundary value problem of system ([1]) with the initial data 

(23) {p,u,E){0,x) = ipo,Uo,Eo)ix), 
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and the boundary condition 



(24) (p,n,E)(t,0) = (pi,-,Ei), p{t,L)=pr 

Pi 

where pi, Ei and pr are the same as in ( |T9|) . 

We assume that the initial data are of the form 



(25) {pq,uq){x) 



and 



{Po-,uq^){x), if < X < xo, 
(Po+,Mo+)(a;), if xo < X < L, 



(26) Eo{x) =Ei+ [ (pois) - b{s))ds 

Jo 

which is a small perturbation of (p, u, E) in the sense that 



(27) 



ko - ^ol + ||(Po+,^^o+) - (p+,m+)||h'=+2([£o,l]) 

+ ||(P0-,%-) - (P-,^^-)llHfe+2([0,£o]) < 



for some small e > 0, and some integer k > 15, where xq = min{xo,xo} and xq = 
max{xo, Xq}. Moreover, (po, Uq, Eq) is assumed to satisfy the Rankine-Hogoniot conditions 
as X = •^'O? 

((p(Po+) + P0+M0+ - (p(po-) + Po-Mo_)) ■ (po+ - Po-)(£o) 

(28) 

=(po+Mo+ - Po-Mo-)^(a;o)- 

Before stating our dynamical stability result, we give the definition for piecewise smooth 
entropy solution for Euler-Poisson equations. 

Definition 1. // (p_, £"_) and (p+, u+, £■+) are smooth solutions of Euler-Poisson 
equations (Q]) in the regions {{t,x)\t > 0, < x < s{t)} and {{t,x)\t > 0, s(t) < x < L}, 
respectively. Then 



(29) (p,M,E) 



(p_,u_,E_), if 0<x<s(t), 
(p+,u+,E+), if s{t) <x < L 



is said to be a piecewise smooth entropy solution of (QP at x = s{t) if (p, u, E) satisfies 
Rankine-Hugoniot conditions 

' {p{p) + pu^){t, s{t)+) - (pip) + pu')it, sit)-) 

= {pu{t,sit)+) - puit,s{t)-))s{t) 

pu{t, s{t)+) - puit, s{t)-) = {p{t, s{t)+) - p{t, s{t)-))s{t), 

Eisit) + ,t) = E{s{t)-,t), 

and Lax geometric entropy condition 



(30) 



(31) 



(u - ^Ap)){t, s{t)-) > sit) >iu- V?G^)(t, sit)+), 
{u+^/^){t,s{t)+)>s{t). 

The dynamical stability theorem in this paper is as follows. 



Theorem 2. Let {p,u,E) be a steady transonic shock solution to system ^\) satisfying 
p8\} . (23); ^2^), ^2^, and ^Mj- Moreover, there exists a 5 > (5 depends on (p, u, E) ) 
such that 

(32) ^_(a;o) = ^+(xo) > -5. 

Then there exists an Eq > such that for any e < Eq, if the initial data {pq,Uq, Eq) 



satisfy /\2o]l . (27), f2o\) and the k + 2-th order compatibility conditions at x = 0, x = Xq 
and X = L, then the initial boundary value problem (Q]), l^[2^ and l^{24\ ) admits a unique 
piecewise smooth entropy solution {p,u, E){x,t) for {t,x) G [0, oo) x [0,L] containing a 
single transonic shock x = s(t) (0 < s(t) < L) with s(0) = xq. 
Furthermore, there exist Tq > and A > such that 



(33) 



{p-,u^, E_)(t, x) = {p^,u^, E_){x), for < x < s(t), 



for t > Tq and 
(34) 

\\iP+^u+,){-,t) - (p+,'"+)(-)llvyfc-7.°°(s(t),L) + \\E+{-,t) - £'+(•) II vF'=-e.°°(s(t),L) < Cee'^'' 



fc-6 



(35) 



J2\dr{s{t)-xo)\<Cee 



xt 



m=0 



for t > 0, where {p±,u±, E±) are the solutions of the Euler-Poisson equations in the 
associated regions. 



8 



The condition f p2|) is used to prove the exponential dynamical stability of the steady 
transonic shock solutions. When this condition is violated, we have the following linear 
instability results for some special cases. 

Theorem 3. There exist L > and a linearly unstable transonic shock solution (p, u, E) 
satisfying and ^ and 

(36) E^{xo) = E+{xo) < -C 

for some positive constant C . 

Begin with several remarks concerning Theorem [1] , Theorem [21 and Theorem [3l 

Remark 2. In both TheoremUl and Theorem{^ the results are also true if we impose 
small perturbations for the boundary conditions 

Remark 3. It follows from the results in [25] and TheoremUl that the background tran- 
sonic shock solution in Theorem [H does exist. Moreover, we do not assume that b{x) is 
a small perturbation of a constant in both Theorem U\ and Theorem [H where it may have 
large variation. 

Remark 4. In [36\, the local-in-time stability of transonic shock solutions for the Cauchy 
problem of (QP is considered by assuming the existence of steady transonic shocks. Here, 
we prove the global-in -time exponential stability for the initial boundary value problem. 

Remark 5. The compatibility conditions for the initial boundary value problems for hy- 
perbolic equations were discussed in detail in [331 ISSl [29] . 

Remark 6. In Theorem IE, the regularity assumption is not optimal. By adapting the 
method in [29], less regularity assumption than that in (21) will be enough . However, our 



proof only involves the elementary weighted energy estimates rather than paradifferential 
calculus. 

As far as the Euler-Poisson equations are concerned, there have been many studies 
on the large time behavior of solutions, see [131 El] and references therein. However, all 
these studies are for the Euler-Poisson equations with relaxations, where the relaxation 
has very strong dissipative effects. 
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It is interesting to compare these results with the transonic solutions of a quasi-one- 
dimensional gas flow through a nozzle. The time-dependent equations for the quasi-one- 
dimensional isentropic nozzle flow are 



where p, u and p denote respectively the density, velocity and pressure, A{x) is the cross- 
sectional area of the nozzle. In [7], steady state solutions for fl37|) containing transonic 
shocks were constructed for the boundary value problem in the interval [0, 1] with the 
boundary conditions (p, u){0) = {pi, ui) and (p, u){l) = {pr, Ur), where (pz, ui) is supersonic 
and {pr,Ur) is subsonic and satisfles piui = prUr- The general wave patterns for flows in 
quasi-one-dimensional nozzles were studied in [21]. The stability of transonic shocks for 
system fl37|) was studied in [22j. It was shown that the transonic shock solution is stable 
if A' > in [22] via a modifled Glimm scheme by introducing some steady states in the 
building blocks. For the studies on the solutions of general hyperbolic conservation laws 
with moving source, see [2T1 [231 1201 [HI [12] and references therein. For piecewise smooth 
initial data, when the nozzle is a straight slowly increasing sectional nozzle, Xin and Yin 
obtained the dynamical stability of steady weak transonic shock solutions, see [lO]. 

Remark 7. The idea for the proof of Theorem [H can also be applied for the stability of 
transonic shock solutions with arbitrary amplitude for quasi- one- dimensional nozzle flows 
with rapidly increasing nozzle walls, see |35j . 

It would be interesting to extend the results in this paper to the multi-dimensional 
case, as those for the gas dynamics, see [21 13S1 [13 [IS] and references therein. An effort in 
this direction was made in [lOj for a viscous approximation of transonic solutions in 2-d 
case for the equations of semiconductors. However, zero viscosity limit in [TU] remains an 
open problem. 

The proof of persistence of transonic shocks under perturbation of b{x) depends on a 
monotonic dependence of the shock location as a function of downstream density proved in 
1121 The a priori estimates for supersonic and subsonic solutions and existence of transonic 
solutions are presented in ^ The proof that for initial data close to steady transonic 
shock the solution decays exponentially to the shock depends on the fact that 

• perturbations to the left as swept upstream by the supersonic hypothesis, and 
importantly that 



(37) 
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• perturbations to the right decay because of absorption at the shock. 

The proof of the latter property is not straightforward. In addition to the usual technical 
difficulties from the quasilinear structure there is fundamental difficulty that the problem 
involves a free boundary (shock) on the left of the subsonic region. The key is to prove 
decay for the linearized problem. After a nontrivial transformation that the linearized 
problem resembles a Klein-Gordon equation ( §4.ip . For that equation we construct a 
nontrivial energy functional which is decreasing due to dissipative boundary conditions 
( §4.2p . Decay to the linearized problem is proved by employing an argument of Rauch- 
Taylor [31] . Then the quasilinear technical argument comes in §4.31 The linear instability 
of the transonic shocks is investigated in ^when the condition fl32l) is violated. Finally, 
we have an appendix for the Ranch- Taylor type estimate |Mj for the linearized problem 
for unsteady Euler-Poisson equations. 

2. Monotone Relation between the Shock Position and the Downstream 



In this section, we derive a monotone relation between the shock position and the 
downstream density (pressure) for the steady transonic shock solutions for the Euler- 
Poisson equations. 

For any supersonic state (p, E) satisfying p < ps, one can connect it to a unique subsonic 
state {s{p),E) via a transonic shock, where s{p) is determined by the entropy condition 
and the Rankine-Hugoniot condition: 



Density 



(38) 



s(p) > Ps and p{5{p)) + 



p{p) + 



P 



This yields 



(39) 



ds 
dp 



(P) 



P\p) 



p'{5{p)) 



Note that smooth solutions of ([8]) satisfy 



(40) 
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Hence the supersonic solution {p{x), E{x)) of (j8]) satisfies 



(41) ^K^^^O _ PE 



The monotone relation between the shock position and the downstream density is given 
by the following lemma. 

Lemma 4. Let {p^^\ E^^^) and (p^'^\ E^'^^) be two transonic shock solutions of and 
(p^^\ E^'^^){i = 1,2) are defined as follows 



(42) {p^'\E^'^] 



iPslp.Ei%), for 0<x< Xi, 
(p?1^12J, for x,<x<L, 



where 

(43) P%<Ps<P^l for ^ = 1,2. 
Moreover, they satisfy the same upstream boundary conditions, 

(44) p(i)(0) = p(2)(o) = p,, i?(i)(0) = i?(2)(o) = El. 
Ifb< ps, xi < X2 and E^upixi) > 0, then 

(45) p«(L)>p(2)(L). 

Proof: Since (pilp, Eiup) and (piup, E^up) satisfy the same ODE system and initial values, 

(46) {p% E%) = {p%, Ei%) for x G [0, x,]. 
For Xi < X < X2, define a function E^ as follows 



(47) 



^ = ^ipfJp) - b, for xi<x< X2, 



E„(xi) = Eiiiix,) = Eiilix,) = EfXixi). 
Note that pstp < ps < 5{pfup), therefore, 
(48) Ef2p{x) < E^{x) for x G (xi,X2]. 
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Since Ea{xi) > and b < p^, one has Ea{x) > for x G {xi,X2]. For x G {xi,X2), the 

/ON 

equation (1411) for ps„p becomes 



(49) 

Thus 

(50) 



rf5(p, 



(2) N 
sup ) 



PsupEsup 



dx 



J2 



5(pili)(a:i) = pitl(a:i). 



ds{pfup) 



< 



PsupEf^ (x) 



P'(5(pia)) - 



It follows from the fact Ea{x) > that, on [xi,a;2], {5{psup), Ea) satisfies 



(2) 



(51) 



ds{p. 



(2) X 
sup ) 



s{piJp)Ea 



dx 



< 



J2 



\Psup) 



dEa 
dx 



Note that (Psti' ^lub) satisfies the equations 



dp; 



(52) 



sub 



dx 



P'ipiub) (Jl) N2 



dE 



(1) 

sub 



dx 



The comparison principles for ODE systems ([30j), yields that 

(53) 4P%)ix2) < Pi'iix2), K(X2) < i?lti(x2). 

Since Eiup{x2) < Ea{x2), we have 

(54) Efi{x2) = Eilix2) < E^{X2) < Ei'J,ix2). 

Note that (pii^, -E^il,) and (pi^6, -E^^),) solve the same ODE system on [x2,i^], by the 
comparison principle for ODEs again, one has 

(55) piti(L) > pfj,(L) and E« (L) > eI'J,{L). 



This finishes the proof of the Lemma. 



13 

□ 



(58) 



3. A Priori Estimates and Existence of Steady Transonic Shock 

Solutions 

In this section, we prove a priori estimates for supersonic and subsonic flows via the 
multiplier method, which yield the existence of supersonic, subsonic, and transonic shock 
solutions. 

A smooth solution of ([8]), p, satisfies the following second order equation: 

(56) (/(p)p')' = p-6, 
where 

(57) /(p) = i (p'(p) - ^) . 

Suppose that (po, Eq) is a supersonic or subsonic solution of steady Euler-Possion equa- 
tions with background charge bo and with initial data {pj,Ei), namely, 

f{p)p = E, 
E' = p-bo 

p{a) = Pi, E{a) = Ej. 

The following lemma gives some stability estimates for both the supersonic and the 
subsonic solutions of (IS^ . which are small perturbations of the solutions of (155]) . 

Lemma 5. For any interval [a, I] C [0, L], let (po, Eq) be a supersonic or subsonic solution 
of the problem ^58\) . Then there exists e > such that if 

(59) \\b{x) - feollcoM + \pi\ + \Ei\ < e, 

then there exists a unique supersonic or subsonic solution {p,E){x) for x G [a, Z] of the 
problem ^ with initial conditions 

(60) p{a)=pi + pi, E{a) = Ei + Ei, 
Moreover, (p, E) satisfies the following estimate 

(61) ||p-po||ciM <Ce"\ 
for some constants C > and a > 0. 
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Proof: We only give the proof for the case when po is subsonic on [a, /] (the case when 
when Po is supersonic is completely similar). When po is subsonic on [a, I], there are 
constants cq > and Ci > such that 

(62) f{po){x) > Co, for x G [a,l], 
and 

(63) Ps < Po{x) < Ci, |p'o(x)| < Ci, for x E [a,l]. 
We may assume a priorily that 

(64) /(p)(x) >Co/2, forxG[a,/], 
and 

(65) ps < p{x) <2Ci, \p'{x)\ <2Ci, ioi xe[a,l]. 

Once we obtain the estimate fl6T]) . under this a priori assumption, the lemma can be 
proved by using the standard local existence theory of ODE systems and continuation 
argument. 
Let us denote 

(66) p = p - Po, b = b-bo. 
Then, 

(67) ifip)p' + Fipo,p)p'oPy-p = -b, 
where 

(68) F{po,p)= f f\p, + ep)de. 

Jo 

Choosing a multiplier K{x) = e~^^^~'"'^ for some constant A > and multiplying both 
sides of fl67|) by K{p' + p), after some straightforward calculations, one has 

f -hK{p' + p)dx = f{{f{p)p' + F(po, p)p'op)' - p}K{p' + p)dx 

J a J a 



(69) 



Kfip) + K[{f{p)y + F{Porp)p'oWY + Kp{f{p)p')'dx 



{K{F{po, p)p',)' + KF{po, p)p', - K) pp'dx 



+ / {K{F{porp)p'o)' -K)p'dx 



Integration by parts and the fact that K'{x) = —XK{x) and K"{x) = A^-ft'(x) yield 
-I 



(70) 



—bK{x){p + p)dx 

I^K{x) ((^ - l) fip) + lifip))' + F{Porp)p'o^ {P'fdx 
+ j'^ K{x) - + {Fp',y - 1)) p'dx 

+ I K{x){iF{po,p)p',y + F{po,p)p'o-l)pp'dx 

J a 



+ K{1) ( + fpp' + A/f 



, - (/^ + fpp + A/—, 

x=l Z Z 



In view of ( 16^ and ( 165|) . we can choose A sufficient large so that 
<C f e-^^^'-^^dx + C{\\pi\^ + \Ej\^). 

J a 



(71) 



Here and in the following, C denotes a generic positive constant C. 
Therefore, if ( 159|) holds for some e > 0, then 

(72) f e-'^^'^'^^dx + Ci\\pi\^ + \Ei\^) < Ct\ 

J a 

This, together with ( 17T|) . implies 



f -p" + {p' fdx < Ce^^'-'^^e^ < Ce^^e' 

J a 



(73) 

which in turn implies 

(74) IIpIIcom < Ce^'^e, 
where a = A/2. Using (167|) . we obtain that 

(75) IIp'IIcom < Ce"^e. 



We are now in position to prove give Theorem [H 
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By the assumption on the unperturbed transonic shock solution (p*-*^-*, -E*-^-') for the case 
when b{x) = Bq^x) (x G [0,L]), there exists a constant 6 > satisfying [xq — 6,xo + S] G 
(0, L), such that the initial value problem for the ODE system: 



(76) 



(p(p) + y)x = pE, E^= p- bo 
{p,E)l=o = {puEi) 



admits a unique smooth solution {p^,E^){x) on the interval x G [0,Xo + S] satisfying 
< p^{x) < ps for X G [0, Xq + S] and 

(77) E^{x) > 0, ioT xo - 6 < X < Xq + 6, 

where xq is the shock location for the unperturbed background transonic shock solution 
(p^^\ E^'^^) for the case when b{x) = bo (x E [0,L]). Moreover, by the uniqueness for the 
initial value problems of ODEs, we have 

(78) (p^ E\x) = (p(°), E(°))(x), for x G [0, Xo). 

Now we define two functions: {p{,E[){x) for x E [xq — 6,L] and {p2,E2){x) for x G 
[xq + 6,L] as follows. Let {p\, El){x) be the solution of the initial value problem of the 
following ODE, 



(79) 



(p(p) + 'j)'x = pE, E^ = p- bo, for X e {xq- 5, L] 



(p, E)ixo -S) = (s(p^(xo - 5)), E\xo - 5)) 



where s is the function defined in fl38|) . It follows from ([5]) that the problem fl79|) admits 
a unique smooth subsonic solution (p^,i?[)(x) on the interval x G [xq — 5,L\ satisfying 
p^(x) > Ps and El{x){x) > E^{xo — 5) > for x G (xq — S, L]. 

Let (p'2, E2){x) be the solution of the initial value problem of the following ODE, 



(pip) + 'j)x = pE, E^ = p-bo,ioi xG(xo + 5,L], 
(p, E){xo + S) = (5(/(xo + S)), E\xo + 6)), 



(80) 



where s is the function defined in (]38l) . Again, following from ([5]) that the problem ( IHOl) 
admits a unique smooth subsonic solution (p2,-E'2)(x) on the interval x G [xq + S,L] 
satisfying P2(x) > ps and E2{x){x) > E^{xo + 5) > for x G (xq + 6,L]. 

By the monotonicity of the shock location w.r.t. the exit density (Lemma Hj), we have 

(81) p;{L) < pr < p\{L). 
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Now, we define two transonic solutions for the case when 6 is a small perturbation of 
60 as follows: Let xi = xo — 6 and X2 = xo + 6 and 



(82) (P«,^«)(x) 



{pi,El){x), ,forO<x<a;i, 
iPi^El){x) ,for Xi<x<L, 

for i = 1,2, where (p^, El){x) is the solution of the following problem: 

{p{p) + 'j)x = pE, E^ = p-b{x), 0<x<Xi 
ip,E)U=o = iPi,Ei) 

and (p[, £■[) is the solution for the following problem 

(p(p) + y)x = pE, E^ = p- b{x), for x E (xi, L] 
{p,E){xD = {5{pl{x.-)),El{x{)) 



13) 



U) 



where s the function defined in ( ]38l) . 

It follows from Lemma [5] that, if \\b — &o||co[o,l] = e is sufficiently small, then {pi. El) 
and {pI, El) are well-defined and satisfying: 

sup rplEl)-{pl,E')\<Ce, 
sup m,El)-{p\,W')\<Ce, 

x(i{xi,L\ 

for 2 = 1,2, for some constant C > 0. In particular, 

(85) |(pI(L)-p[(L)|<Ce, 
for i = 1,2. This, together with (IHTi) . implies that 

(86) pliL) < pr < p{{L), 

if e is sufficiently small. This shows that the transonic shock solution of the boundary 
problem ^ and ([2]) admits a unique transonic shock solution (p, E) with a single transonic 
shock located at some Xq G {xi,X2) by a monotonicity argument as follows: For a G 
[a;i,a;2], we define a function (^(a) = p(L) where p is a transonic shock solution of the 
system ([8]) with p(0) = p; and -E(O) = -E; and its shock is at a. By Lemma[5l we know that 
g{a) is a continuous on [xi, X2]. Moreover, Lemma H] implies that g{a) is strictly decreasing 
on [xi,X2]. Finally, the stability estimate. Lemma O yields that xq G [xq — Ce,xo + Ce] 
for some constant C > 0. 
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Proof of the uniqueness of transonic shock solutions. Suppose, besides the transonic 
shock solution we constructed above, there is another transonic shock solution for the 
perturbed boundary value problem ([8]) and ([9]) for the variable b{x), say, {pi,Ei) , which 
contains a single transonic shock at x = Xi. Because of the monotonicity of the shock 
location w.r.t the downstream density, xi must lie out of the region (a;o — 6,xo + 6). Let 
g be the function defined on [0, L] as g{a) = p{L) where p is a transonic shock solution of 
the system (|8]) with p(0) = pi and E{0) = Ei and its shock locates at a for the case when 
b{x) = bo. By Lemma [5l g is a. continuous function, so 

(87) S = {y\y = g{a), a G [0, L]\{xo -6,Xo + S)} 

is a closed set. Since we assume that the transonic shock solution (pi, ^i) to the problem 
(8) and (9) for the case when b{x) = 6o is unique, pr ^ S. Hence, there is an ei > such 
that 

(88) inf \pr — y\ > ei. 

yes 

In particular, 

(89) -^(xi)| > ei. 
On the other hand, by Lemma we have 

(90) \pr - g{S:i)\ = \pr - PiiL)\ = |pi(L) - pi(L)| < ei/2. 

if b{x) satisfies f|T71) for some eo > 0. This contradicts fl89l) . Hence, there is only one 
transonic shock solution for the problem ([8]) and (Q. 

4. Dynamical Stability of Transonic Shock Solutions 

In this section, we investigate the dynamical stability of transonic shock solutions for 
the Euler-Poisson equations ([T]). 

4.1. Formulation of the Problem. Let (p, m, -E) be a steady transonic shock solution 
of the form (ITS!) satisfying (l32l) . Suppose that the initial data (po, uq, Eq) satisfies (l27|l and 
the /c + 2-th order compatibility conditions. It follows from the argument in [12] that there 
exists a piecewise smooth solution containing a single shock x = s{t) (with (s(0) = Xq) 
satisfying the Rankine-Hogoniot conditions and Lax geometric shock condition fl30|) . of 
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the Euler-Poisson equations on [0, T] for some T > 0, which can be written as 



(91) 



{p,u,E){x,t) 



{p^,u^,E^), if < X < s(t), 
{p+,u+, E^), if s(t) < X < L. 



Note that, when t > Tq for some Tq > 0, {p^,u^, E_) will depend only on the boundary 
conditions at x = 0. Moreover, when e is small, by the standard lifespan argument, we 
have Tq < T (cf. [H]). Therefore, 



In the following, for simplicity of the presentation, we may assume Tq = without 
loss of generality. We want to extend the local-in-time solution to all t > 0. In view of 
(!92|) . we only need to obtain uniform estimates in the region x > s{t),t > 0. For this 
purpose, we will formulate an initial boundary value problem in this region. First, the 
Rankine-Hugoniot conditions for (1911) read 



(92) 



(p„, M_, E_) = (p_, u_, E_) for t > Tq. 



(93) 



[pu] = [p]s'{t), [pu^+p\ = [pu]s'{t), 



where [/] = f{s{t)+,t) - f{s{t)-,t), so 



[p + pu^] ■ [p] = [pu]'^. 



More precisely, one has 



p(p+)(t, s{t)) + -^{t, s{t)) - p(p_)(t, s{t)) 




{J+{t,s{t))-J4t,s{t))f, 



where J = pu. In view of ( l92i) . 



(94) 



J{s{t)-,t) = J. 



Hence, 




p+ 




{J+{t,s{t))-Jf. 
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It follows from the Rankine-Hugoniot conditions for fl2Tl) and Taylor expansions that 

p'ip+){s{t)){p+{t, s{t)) - p+is{t))) - ^{s{t)) ■ (p+(t, s{t)) - p+ism 

p+ 

2J - P 

+ —{sit)) ■ (J+(t, s{t)) - Jis{t))) + dMP+) + — )(a;o) ■ {s(t) - Xo) 
p+ p+ 

J- 1 - - 

-d^(p{p-) + —){xo) ■ - a;o) + Ri > {p+{xo) - p-{xo) + R2) 
P- } 

= (J^.(t,.(t))-J>(t))f, 

where 

R, = 0(l)((p+ - p+)2 + (J+ - J)' + (sit) - Xo)') . 

Later on, we always use -Ri to denote those quadratic terms with different 0(1) coefficients. 
It may change from line to line. Thus, by implicit function theorem, we have 

(95) (J+ - J){t, s{t)) = M((p+ - p+)(t, s{t)), s{t) - Xo) 

where s^/i regarded as a function of two variables satisfies 

n , ds^i {p+-p.)E+, , 

M0,0 =0, jr, '—^= -j—^{xo), J- S = - i^o). 

d{p+-p+) ^J/p+ d{s-xo) 2J/P+ 

Substituting fl95l) into the first equation in fl93|l yields 

(96) s'(t) = ^/2(p+-p+,s(t)-xo) 
where £^2 satisfies ^2^(0, 0) = and 

d^2 ^ p'ip+)-P/pl ^^^^ d^2 ^ -i±{xo) 

d{p+-p+) 2m+(p+-p_) ' d{s{t)-xo) 2u+ 
It follows from the third equation in ([1]) that 

E+{x,t) = Ei+ / (p_-b)iy)dy+ / {p^-b)iy)dy, 

Jo J s(t) 

for s{t) < X < L. Applying the first equation in ([1]) and the Rankine -Hugoniot conditions 
( 1^ . one has 

dtE+ = -p+u+{t, x) + piui = -J+{t, x) + J. 

Set 

(97) Y = E+{x,t) - E+{x). 
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Then 

(98) Yt = J-J+, Y, = p+-p+. 

Therefore, it follows from the second equation in the Euler-Poisson equations ([1]) that 

- E^)u + (p{p+) + — - p(p+) -^)+ - = 0. 

V P+ P+J 

Thus, 

(99) duY + fp(p+) + — - p(p+ + Y^) - + E+d^Y + p+Y + YY, = 0. 
One can write this equation in the following way. With ^ = (^o,^i) = and 



for ij = 0, 1, 

(100) Yl a^A^,Yt,Y^)^,,Y + J2bi{^^Yt,Y^)d,Y + c{x,Yt,Y^)Y = 0, 

0<i,j<l 0<i<l 

where ciij, hi and c are smooth functions of their arguments, and satisfy 

^^Y = ai,(x,0,0)9i,F+ ^ 6i(x,0,0)9iy + c(x,0,0)y 

0<j,i<l 0<i<l 

101 - --- 

P 2 J - 

= duY - dM{P+) - —)d:,Y) + d^{—dtY) + E+d^Y + p+Y. 

P+ P+ 



Furthermore, the Rankine-Hugoniot conditions fl95|) and fl96|) can be written as 

(102) rt = -M(v;,s(t)-xo), 

and 

(103) s' = £/2{Y^,S-Xq), 

respectively. Moreover, direct computation yields 

Y{s{t),t) = E^{s{t),t) - E^{s{t)) 
= E.{s{t)) - E^{s{t)) 

= E_{xo) + E_{s{t)) - E_{xo) - E+{xo) - E+{s{t)) + E+{xo) 
= d:rE^{xo) ■ (s(t) - xo) - d^E+{xo) ■ (s(t) - Xo) + Rl. 
Using the third equation in the Euler-Poisson equations ([T]), one has 

(104) s{t)-xo = ^3{Y{t,s{t))), 
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where =2^3(0) =0 and 



1 



dY p_(xo) - p+(xo)' 
Combining f ll02p and f ll04p together yields 

(105) dtY = £/4{Y^, Y), at X = s{t), 

where 



(106) ^.(0,0) = 0, ^^^1^^, , ^^^i^^j. 
Note that on the right boundary, x = L, Y satisfies 

(107) d^Y = 0, at X = L. 

Our goal is to derive uniform estimates for Y and s which satisfy (llOOp . (11041) (llOSp and 

(unzD- 

Introduce the transformation 

(108) t = t, X = {L - xo)j—^ + xo, a{t) = s{t) - xo, 
to transform the problem to the fixed domain [a;o,i^]. Let 

(109) q,ix,a) = ^ ^ ~ - , q2{a) ^ ~ 



L — xo — cr{t) L — xo — cr(t) 

It is easy to verify that 

d d , ~ Odd 
dt di ^ '^^dx' dx 



Qt2 ' '''dx^ ''dxdi ^ V L-XQ-a{i)J dx' 

dxdt '^^ \dxdi L — xo — a(t)dx ^^"^ 55;^/' ^^55;^ 
So (!99|) becomes 



2or'(t)gi%F + (gi(T'(t))29^£F + p+F + ^+52^^^ + 52!^^ 



L — xq — a(t) 
--a"{t)q^dsY. 



rrQidxY 
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With the help of straightforward calculations, the equation fll04p becomes 
(110) a = £/3{Y{t,x = xo)). 

The equation for the shock front, (llU3p . becomes 

(111) 

Applying f lllOp to represent the quadratic terms for a in terms of F, we have, ai x = xo 
(112) 

where ^2 satisfies 

<C{Yi + Y'). 



i = ^2(g2(cT)y^, a{t)). 
at 



^ + ^{xo)a = ^,{Y„Y), 
at 2u, 



%{Y„Y)+ 'y^K 5 ixo)Y, 



2(p+ -p_)u+ 

Clearly, it follows from f lllOp and f lll2p that one can represent cr and cr' in terms of Y and 
its derivatives at x = a;o. Thus, after manipulating (11051) with (IllOp and (I112p . we have 

(113) Yi = ^i{Y,,Y), atx = xo. 

Or, by the implicit function theorem again, equivalently, we have 

(114) Y^ = ^,{Yi,Y), atx = Xo, 
where ^3 satisfies 



(115) 



%{Yi,Y)- 



2u, 



■(^o)n 



;ii6) 



In the following, we drop ~ in x and t for simplicity of notation. 
The problem takes the following compact form 

££{x, F, a)Y = a"{t)q,d,Y, (t, x) G [0, 00) x [xq, L], 

9,y = di{Yt, Y)Yt + e.iYt, F)F, at x = xq, 

dxY = 0, at X = L, 

a{t) = £/s{Y{t,Xo)), 

where, by using ^0 and ^1 to denote t and x, respectively, 

1 1 
^(x, Y, a)Z = y aij{x, Y, VY, a, (j')dijZ + V ^^(x, F, VF, a, a')diZ 
(117) ,^0 



i=0 



+ ^/(x,r,VF,a,a')Z, 



24 

with 



(118) d,{Y„Y) 



Furthermore, one has 



dY, 



m,eY)d9, e,{Y„Y) 



dY 



{9Yt, 9Y)de. 



(119) 



and 



(120) 



^(x,0,0)Z = ^qZ, 



aoo(a;, Y, VF, cr, a') = 1, aoi(x, 0, 0, 0, 0) = aio(x, 0, 0, 0, 0) = — = m+, 
au{x, 0,0,0,0) = - (^p'{p+) - 

bo{x, 0, 0, 0, 0) = d,, (^^^ , hix, 0, 0, 0, 0) = -d, (p'ip+) - ^ 1 + 



gix,0,0,0,0) = p+, di{0,0) 



2n_ 



c2(p+) - 



(xo), ei(0,0) 



^4 



C2(p+) - Ul 



(Xo). 



4.2. Linear Estimate. In this subsection, we study the hnearized problem. 



Theorem 6. Assume that satisfies nS^I . Let Y be a smooth solution of the linearized 
problem 



(121) 



' if (x, 0, 0)F = 0, xo < X < L, t > 
d.Y 



u-^^^ -2 i^o)dxY + — 2-(xo)F, atx = xq. 



dxY = 0, at X = L, 
^ Y{0,x) = hi{x), Yt{0,x) = h2{x), xq < x < L. 

Then there exist G (0, 1) and T > such that 



(122) 



^k{Y,t + T) <ao0k{Y,t) fort>0. 



where 0k is defined as follows 

k 

(123) 0k{Y,t) = Y,^m{y,t) 



m=0 
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with 



(124) 



^uy,t) = —ixo)idrY)\t,xo) 

p+ 



1 

p+ 



{dr'Yf + ( p'(p+) - ) {d^dTYf + p+^rr^ (t, x)dx. 



Remark 8. When Ej^ satisfies ^3^ . it follows from the Sobolev inequality that 
(125) ^m{y,t) > C r {{dT^^Yf + {d^dTYf + dTY^] {t,x)dx. 

Jxo 

This is the key reason that we can handle the case that E^{xq) is a negative number with 
small magnitude. 

Proof: The proof of the theorem has four steps. 

Step 1: Energy estimates. Multiplying the first equation in fll2ip with -^^^-^ dtY on both 
sides and integrating by parts, we have 



=-dt 
2 



+ 



L 



-d. {-)- + -d. ] ] {d,Yfdx 
\P+J p+ p+ \p+. 



(126) 



+ 



J 



+ ^{dtYY 



XO P+ 



P\P+) - ir ) dtYd^Y 



Noting that J is a constant, and using the equation ([2]), we have I2 = h = 0. Apply the 
boundary conditions f ll2ip to get 



h ~{dtYf 



J 

--L Pi 



idtYf 



x=xo 



(127) + - (p'ip^) - ^) d,Y ( '^''\,,-, d^Y 

P+\ P+J \P {p+) - J^IP+ 



J 



{dtYf 



J 

+ — 

x=L p\_ 



{dtYf 



x=xo 



a 1 14^1- 

2p+(a;o) 



p'(p+) - jyp^ 



x=xo 



-Y 



X=Xq 



26 

Thus, 

(128) ^o{Y,t) + D^{Y,t) = ^^{Y,Q), 
where 

(129) Do{Y,t) = 2([ L{dtY)\s,L)ds+ f ^{dtY)\s,xo)ds 



Step 2: Rauch- Taylor type estimates. Using the boundary condition at x = xq and the 
fact that ^^2u+^^ (••^o) — ^ some constant C > 0, it is easy to see that 

(130) Do{Y,t)>Ci [ {Y^ + Y^){s,xo)ds-C2 I Y\s,xo)ds, 

Jo Jo 

for some positive constants Ci and C2, independent of t. Therefore, 

(131) MY,t) + Ci [ {Y,^ + Y^){s,xo)ds < ^o{Y,0) + C2 [ Y\s,x^)ds 



Following from the argument in [3l] and the details in Section [6l there exists a T > and 
5 G (0, r/4) such that 

(132) r{Y^^ + Y^){t,x^)dt> ['''\o{Y,s)ds-Cs f Y\t,xo)dt. 

Jo J^-s Jo 

Since ipQ{Y,t) is decreasing with respect to t, in view of fll28p . we have 



(133) 



(r,2 + Y^){t, xo)dt >6My, \ + 5)-C^ I Y\t, xo)dt 



>6MY,T)-C3 [ Y\t,xo)dt. 
Jo 

Combining f ll28p . f ll30p and f ll33p together, we have 

(134) (1 + C4)^o(V', T) < ^0(1^, 0) + C, f Y^{t, xo)dt, 

Jo 

for some positive constants C4 and C5, independent of t. 

Step 3: Spectrum of the solution operator. Let us define a new norm || • ||x for the 
function h = {hi,h2) e x L'^{[xo,L]), 



(135) \\hf^ = ^(xo)\h\'ixo)+ / —{\h2\' + ip'ip+)-^)\h[\' + p+\h\'\ix)dx. 

p+ Jxo p+ [ p1 

The associated complex Hilbert space will be denoted by (X, || ■ ||x)- Define the solution 
operator 6*4 : X i—)- X as 

(136) Stih) = iYit,-),Ytit,-)) 
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where Y is the solution of the problem (11211) with the initial data h = {hi, h2). By (11281) . 
we can see that St is bounded and satisfies 

(137) ll^ill < 1. 



It follows from the spectrum radius theorem [TB] that the radius of the spectrum of St is less 
than or equal to 1, i.e., |cr(S'i)| < 1. Furthermore, we can define a map i^T : X i— )■ lv^([0, T]) 
as 

(138) K{h) = Y{t,xo). 



Thus the estimate (11341) can be written as 

(139) {l + C,)\\ST{h)U< ||/i||x + C7||ir/i|U2([o,T]), 

for some positive constants Cq and C7. Note that for the initial data /i G X, there exists a 
solution Y G if^([0,T] x [xq,/^]). It follows from the trace theorem of space, [8], that 
K is compact. Applying the Lemma on page 81 in [32j, we have the following proposition: 

Proposition 7. There are only finite generalized eigenvalues for the operator St in the 
annulus {jipc^ < \A — 1} complex plane, each of these eigenvalues has the finite 

multiplicity. 

Step 4: Refine estimate for the spectrum of the solution operator. 
Proposition 8. There is no generalized eigenvalues of St on the circle \z\ = 1. 
Proof of the proposition [8t Suppose not, there exist a; G M and G X such that 

(140) {St - e''^I)V = 0. 

Note that the identity (11281) still holds in the complex setting if we replace the square 
terms in ipo and Dq by the square of modulus. Thus 

(141) ifo{y,0) - My,nT) = nDo{Y,T). 

Since both ipQ{Y,nT) and ipo{Y,0) are both positive and finite, it yields that 

(142) Do{Y,T) = 0. 
Therefore 

(143) ipo{Y, t) = ipo{Y, 0), for all t. 
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Let 



(144) Y = ker{ST - e''^!). 

Note that the coefficients in the problem (11211) do not depend on t, so 

{St - e'''I)St = StiSr - e'^I) 

In particular, S(t)Y C V for any t, so V is invariant with respect to St- Therefore, Stly is 
a semigroup on a finite dimensional subspace. This yields that 

(145) ^(t)|v = e*^ 

for some matrix A. Choose an eigenvector w of A such that 

(146) Aw = Xw. 
So 

(147) Stw = e^'w. 

Note that if St has an eigenvalue on the unit circle, then there is a solution with A pure 
imaginary. If StW = e*^*w, then it follows from fll42p that w{xo) = 0. Furthermore, by 
the boundary condition at Xq, we have w'{xo) = 0. Note that w is an eigenfunction, so 
w = 0. This finishes the proof of the proposition. □ 
Combining Proposition [7] and Proposition [8] together, one has 

(148) cr(^r) C {|z| < for < tto < 1, 

where alSx) is the spectrum of St- Since the coefficients for both the equation and the 
boundary conditions are independent of t, applying this estimate to the time derivatives 
of the solutions implies that 

(149) <^„(F,T) <ao<^„(y,0). 

Summing these estimates together, we have fll22p . 

This finishes the proof of the theorem. □ 



The consequence of the theorem is that the solution of the problem (I12ip decays expo- 
nentially. 
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Corollary 9. There exist constants Aq > and C > such that for any solution Y of 
the problem (QSp, one has 

(150) 0k{Y,t)<Ce-^''0k{Y,O) 
and 

POO j^ + l 

(151) / e^J2^\dlY\\t,xo) + \dlY\'it,L))dt<CMy,0)- 

•^'O 1=1 

Proof: Noting that 0k{y, t) is decreasing in t, for any t G [nT, {n + 1)T), n G N, one has 

0kiY,t) <0k{Y,nT) < a'^0k{Y,O) 

<< <^fc(y,o) = e-"<%-V(>^,o), 

where we have chosen Aq = —^f^. 
Note that 

21+171 fe+i 

XQt 



e^$^(|5iV'p(t,Xo) + |aiFp(t,L))cit 

2^+171 fc+1 



(153) <e^«2»-^^ /" V(|9^Fp(t,Xo) + \dlY\\t,L))dt 

J2'T 

=e^'^'''^{0k{Y, TT) - (^(F, T+^T)) 

<e^«2'-^ V^e-^^'^^,(F, 0) = ao^e-^^'-^^^,(y,0), 

where we use the energy equality fll28p . Thus 

/ e^{\d^,Y\\t,x,) + \d[Y\\t,L))dt 
Jo 

(154) =E / e^5^(|9irp(t,Xo) + |9^Fp(t,L))rft 

00 

<5^«o-V^°^2"">fc(F,o) < c^fc(y,o). 

This finishes the proof of the corollary. □ 
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4.3. Uniform A Priori Estimates. The existence of local-in-time solutions for the 
problem (HD, 1^, 1^, (EO]), (EI]), and ([22]) can be achieved via the argument in [19]. In 
order to get global existence of the nonlinear problem ([1]) , (ITH!) , (fT9l) , (120]) , ([2T]) , and ([22]) , 
we need only to derive global a priori estimate for the problem 

' ^(x, Y, a)Y = a"{t)qi{x, cr)d:,Y, xq < x < L, t > 0, 

d^Y = diiYt,Y)Yt + eiiYt,Y), at x = xq, 

(155) <! (x{t) = j2/^{Y{t,xo)), 

dxY = 0, at X = L, 

Y{0,x) = hi{x), Yt{0,x) = /i2(x), xo < X < L. 

For t > T in Theorem [6] and an integer k > 15, we introduce the following notations 
(156) 



|||(r,a)||||= sup Yl (E^^II^^^^ 











dV^ 



-|- sup 

0<T<t 



\\d\dr'y{r,-)\\L^i\.oM)+ E ll^*^^'"'^(^r)l|L^([.o,L]) 



0<Z<m, 
yo<m<A: 



0<Kfc 



+ sup ||9f+iF(r, 

0<T<t 



0<l<m, 
0<m<fc+l 



0<i<m, 
0<m<fc+l 



df' 



L2[0,i] 



0<m<fc+l 

where A is defined in (I223p . 

In this subsection, when there are no specific indications, we always assume ajj, 6j, and 
g are functions of (x, Y, VF, a, d), (ii and ei are functions of (F, y(). Furthermore, for any 
/ e N and given Y and cr such that || || (F, o") || || < oo, we define 

(157) ^i{Z, t- F, a) = <l;_i(Z, t; F, a) + $o(9^Z - gi(x, ^)>^x^, ^), 



and 
(158) 



^i{Z,t-Y,a) = Y,'^n^{Z,t-Y,a), 



m=0 
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where 
(159) 



By fll20p . it is easy to see that 

(160) ^miZ,t;Y,a)it)>C f\dT^^Zf + {d^dTZf + {dTZf){t,x)dx, 

for some constant C > independent of t if || ||(y, cr)|| || < e for some e > 0. Also, 

(161) <l>fe(Z,t;0,0) = ^fc(Z,t). 

Proposition 10. Assume that satisfies [3^) . There exists an eo > 0, for any e G 
(0, Co), if {Y, a) is a smooth solution of the problem /1155\) with {hi, h2) satisfying 

(162) |o-o| + ||/ii||^fc+2 + ||/i2||i^fc+i < < Cq 
and II II (F, cr)|| || < e for t > T, then 

(163) ||ll(>^,^)ll||<|- 

Proof: We divide the proof into four steps. 

Step 1: Lower order energy estimate. Taking the m-th (0 < m < /c — 1) order derivative 
for the equation (11551) with respect to t, then 

(164) if (x, y, a)drY = ^„(x, Y, a) + #^(x, Y, a), 
where 

(165) ^mix, Y,a)= J2 dla.A,dr-'Y - ^ dibAdr'Y - dlgd^'Y 

l<l<m I i,j=0 i=0 



and 

(166) #^(x,r,a)= J2 CL^dr'iqiix,a)Y^) 

0<l<m 

where and in the following is the binomial coefficient, 

ml 



if m > /, 



cL = l Hm-iy. 

0, if m < I. 
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Multiplying the both sides by d^^^Y/p^{x) and integrating on ^2 =: [0,t] x [a^o,-^'] for 
t e [0,T], noticing that aoo = 1, we get 



/ 



/ / ^ix, Z, * , , drdx 

JJn P+{x) 

[t, x)ax 



(2p+(x)) 

bo o / «01 \ \ /Qm+lv-\2, 



In \P+{x) P+ 



(167) +yX— ^(^'-) - 2 







Jo V / 



-I 



-(0, x)dx 



1=1 

In view of the fact 



&o(x,0,0,0,0) _ ^ ^ aoi(3;,0,0,0,0) 

1 ^aJ^^VaJ^^ 



P+{x) \P+J \ P+ J 

=0 



and 



h{x,0,0,0,0)__^ f an{x,0,0,0,0) 



P+ix) 



1 



P+{x) 



P+ix) 



<9^ 



-(pXp+) - JVp' 
p+ 



we conclude 

\J2 + M <c 

Note that 

(168) 

thus 



2 / I 



/=0 \ i=0 



m+1 1 



m+1 1 



1=1 i=0 



1=1 i=0 



\J2 + Jsl <C sup 
0<r<t 



At 

e«4 



2 / I 



1=0 \ i=0 



X 



(169) 



/' 

Jo 



e el^^(Y,T;Y,a)dT 



<C \\(Y, 



»ll|| / 

^0 



e el^^{Y,T;Y,a)dT 



Similarly, 
(170) 



IJ4I < C^||||(i",<T)|||| f e-^^m{Y,T;Y,a)dr. 
Jo 

It follows from the boundary condition at a; = L, Yx{t, L) = 0, that 

(171) J,= ^{dr^'Y)\T,L)dt>0, 

Jo P+ 



if II II (r, a) II II is sufficiently small. 

Differentiating the boundary condition 



(172) 



d,Y^di{Yt,Y)Yt + e,{Y,,Y)Y 



m times with respect to t, one has 

(173) d'^Y^ = did'^+^Y + e^d^Y + at x = xq, 

where satisfies 



(174) 
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Therefore, at x = Xq, 



(175) 



So 



iandi + aoi){dr^'Yy-dt{anei 



+ dtianCi) 



(dry? 



- dt (aiici) 



p+ 



(t,Xo) 



(r, Xo)c?r + 



P+ 



Xo)dr 



(176) 



> 



(ar+'r)^(r,a;o)rfr 



* — (aiifii + aoi) ^^m+i-i^\2 



aiici- 



m+l 

^1111(^.^)1111 / J2(9lYriT,Xo)dT. 



Summing up all the estimates fll69p . fll70p . (I17ip . and (11760 together yields 

Sf{x, Y, a)d"'Y t , , dtdx 
p+{x) 



>$„(r, t- r, a) - $™(r, O; r, a) + ^„(y, t; y, a) 

e-^<l^(r, r; F, a)rfr + #„(r, t; y, a) 



(177) 

where &m is defined as follows 

(178) ^™ = 5^^K^,^;>^,^) 

1=0 

and 



(179) ^K^,i;>^,^)= 

Jo 



-{audi + gpi) 



(a^^Z)^(r, a:o)rfr + ^(a^^Z)2(r, L)dr. 
Jo P+ 



In particular, ^Q{Z,t]0,0) = Dq^Z^I) which is defined in (11291) . 
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Note that 



(180) 



p+ 



l<l<m, ' l<l<m, 

0<i,j<l 0<i<l 

I. 



+ // E ci^dr'Y&r'Y 



1<1<6, 
0<ij<l 



+ fiJLc'm^^dr'Ydr'Y 



KK6 



7<l<m 
0<i,j<l 



We estimate each term as follows. The first part of Ki can be estimated as follows 



(181) 



<C 



l</<6, 
0<i,i<l 



KK6 



C. 



I dlai. 



'p4x) 



d\a. 



p+[x) 



L°°([xo,L]) 



Using the equation fll55p . we have 
(182) 

m+l / / m+l ji 

EEll^*^^^^(^'-)lli^[^o,L]<C^U™(>^,T;y,a) + ||||(F^ ^ 

1=1 i=0 \ 1=1 
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Thus 



(183) 



l<i<6, 
0<ij<l 



<C\\\\{Y,a 



Xt 

e 64 



m+1 



^m{Y,r;Y,a) + \\\\(Y,a)\\\\J2 



1=0 



IF 



dr. 



The estimates for the second part of Ki and K2 are quite similar to the estimate above. 
Moreover, for the first part of iCs, we have 



• 7<l<m, 
0<i,j<l 



(184) 



^0 7^;^, 



7<l<m 



d[aij 



<ciiii(>^,^)iiii r 

^0 



1(9- • (9™"' y 1 1 



L2([xo,L]) 

^m{y,r;Y,a) + \\\\{Y,a)\\\\J2 



At 
e 64 



rrt+l 





2" 







dr. 



So, we have 
(185) 

grn+ly 



Iln 



n P+ 



dtdx < C\\\\{Y,a 



f 

Jo 



_ At 
e 64 



m+1 



My,r;y,a) + ||||(y,<7)||||5^ 



1=0 





2" 


dt^ 





(ir. 
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Similarly, 



<C 



p+ 



0<K6 



dt^- 



+ c 



7<l<m 



d'+^a 
~di^ 



P+ 

5m-/ „ ^^ t ^^^^ 



"^0 0</<6 



a 



(186) 



7</<m 



JO 7 



dt^+^ 



IL2([xo,L]) 



L2([xo,L]) 



Il°°([xo,L]) 



Ar 

e 64 



m+2 



(it 



+ c||||(r,a)||| 



At 

e 64 



<c^||ll(>^,^)ll|| 

where we have used the estimate 



At 

e 64 



$^(y,r;F,a) + ||||(y,a)|||| J] 

m+2 
m+2 

<l„(F,r;y,a) + J] 



L2([xo,L]) 

dr 



d^a 


2" 


dt^ 





/=0 



d^a 


2" 


dt^ 





d^a 


2" 


dt^ 





dt^^ 



< C 



a 



dt^+^ 



Combining the estimates fll77p . fllSSp and fll86p together, one has 
^Uy, t; Y, a) + t; Y, a) < <I„(F, 0; F, a) 

(187) 









/ e 64 




Jo 


■ ■ ■ ,fc - 1. 





m+2 



l'^(F,r;y,a) + ^ 



z=o 





2" 







dr + ^™(F,t;y,a) 



for m = 0, 1, 

Step 2: The highest order energy estimates. Take the fc-th order derivative for the 
equation fllSSp with respect to t, then 

^{x,Y,a)d^Y = ^k{x,Y,a) 



0<l<k-l 
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In order to handle the term , we rewrite the above equation as 

(189) ^{x, Y, VY, a, &)Y = ^k{x, Y, a) + F, a), 
where 

(190) Y = d^Y -q,{x,a)Y/'''^ 



and 



#(x,r,a)= Ci-—^dt\qiMY^) 



0<l<k-l 



(191) - (^^dt{qi{x,a)Y,) + ^dUqi{x,a)Y,) 



(192) 



1 /^k^ 

- (2aoi<9t<9x + audi + ^ bidi + g) { -j^Qii^^ cr)^^ 

Note that the boundary condition for Y at x = Xq is 

d.Y ={di + ^^Y, + ^Y)d,Y + e,Y 
oYt dYt 

+ + + A^TT 

+ Ecl.okII)^. + a;(|i)y)9f"-'y 

+ dt-\^d,Y)d,Y + dt\^d,Y)Y 

fc-i fe-i 

+ 5^c^9^rfi(F„r)9r+iF + 5^ci9^ei(y„F)9rY. 

1=1 1=0 

The associated boundary condition for F at x = L is 

„ - 9 / L — X d''a\ 

Multiplying the both sides of f ll89p by and integrating on fi, and noting that 

fc+i / fe+i 



(194) $^||5^e'n^,-)||i.[.„,L] < C I I..(y,r;y,a) + ||||(F,a)||||'5^ 

=0 



dt^ 
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one has 



(195) 



%{Y, t- r, o) + %{Y, t- r, a) < $o(r, O; r, a) 

k+l 

^,{Y,T-Z,a) + Y, 

1=0 









/ e 64 




Jo 





2" 


dt^ 





where is defined as follows 



(196) 



dT + ^kiY,t;Y,a) 



d^a 



,(Z, t; F, a) = ^,_i(Z, t; Y, a) + ^o(5f ^ - qi{x, a)F„, t; F, a). 



dt^ 



Adding the estimates (11871) and (11951) yields 

l>fe(y,t;y,a) + #fe(F,t;F,(T) < <lfc(y,0;y,cr 

(197) 







^k{Y,r-Y, 


{fo'~'' 









k+l 



1=0 





2" 


dt^ 





dT + ^k{Y,t-Y,a) 



Step 3: Boundedness of the energy. Differentiating the equation for the shock front 
(198) a{t) = s^^{Y{t,xo)) 

with respect to t, we have 



k+l 



1=0 



(199) E 
Therefore, we have 

(200) 



dt^ 



k+l 

E 

1=0 



k+l 



{r) <C [\a{r)\' + Y.\^lY{r,Xo)\' 



1=0 



a 



dt^ 



< 



1=1 



Combining (I199p and (I200p together, we get 



(201) 



t k+l 

At 
e 64 



=0 



dT<C{ / e-6i<l>o(r,r;F,(T)rfr + 4(>^,t;>^,fx; 



Thus the energy estimate (11971) is equivalent to 

<lfc(r, t- r, a) + %{Y, t- r, a) < <lfc(F, O; r, a) 

(202) 



If II II (F, a) II II < e, then 



(y,r;r,a)+ / e-64$,(r,r;F,a)rfr . 



(203) 



<lfc(F, t; Z, a) + ^fc(r, t; Z, a) < 2<lfc(r, 0; Z, a) < Ce' 
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This yields 



sup 

0<T<t 



0<l<m, 
\Q<m<k 

+ sup ||9f+^F(r, 

0<T<t 



0</<A: 



(204) 



5^ ||9^F(-,Xo)|U^[o,*]+ 5^ \\dldr'Yi;L)h2[o,t] 

l<Kfc+l l<i<fc+l 



+ E 



l<m<fc+l 



L2[0,t] 



<C ( sup <|i/2(r, r; y, a) + ^^/'(r, t; F, a] 

\0<T<t 



<Ce^ < - . 
~ ~ 4 

Step 4: Decay of the lower energy and the shock position. The basic idea to get the 
decay is to control the deviation of the solution Y of the problem f ll55p with the solution 

Y of the problem (I12ip via Y . The contraction of the energy for Y will also yields the 
contraction of the energy for Y . This gives the decay of the lower energy of Y . The decay 
of the shock position is a consequence of the governing equation for the shock front and 
the decay of the lower energy. 

At time r = to, we can choose hi G and h2 G H^^^ such that there exists a solution 

Y e C^-^-\[tQ,oo)]H\[xQ,L\)) of the linear problem f imi) satisfying F(to, •) = and 
Yt{to^-) = h2, and Y satisfies 



(205) 



k-l I 
1=0 1=0 



for some uniform constant C, and 



(206) 
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Note that Y — Y satisfies the equation 



(207) 



ij=0 

1 



j=0 



J2 ^) - ciijix, 0, 0))^^,? + 5^(6,(x, Y, a) - bj{x, 0, 0))d,Y 

i,j=Q j=0 

+ {g{x, Y, a) - g{x, 0, 0))F + a"{t)q,{x, a)dxY, 



and the boundary conditions 



(208) 



d,{Y -Y) = d,{Yt, Y)dt{Y -Y) + ei(Fi, Y){Y - Y) 
+ {di{Yt, Y) - c/i(0, 0))dtY + {e,{Yt, Y) - ei(0, 0))F, at x = xq, 



and 



(209) 



a^(F -Y) = 0, at X = L. 



Taking the m-th (m < fc— 4) order derivatives with respect to t for (12071) and multiplying 



5™+^(y-y) 



(if m < — 4) or 



dtZ 



the both sides of the resulting equation with 
gk-4(j^ - F) - ^^^qi{x,a)Y^ if m = A; - 4), and integrating it on [to, to + T] x [xq, L] 
where T > is the one appearing in Theorem |6l the new terms other than those appearing 
in the energy estimates in Steps 1 and 2 should be induced by the first three terms on 
the right hand side in (120 7p . Here we only give the highest order energy estimates. Let 
us estimate the first one as follows 



pto+T pL _ d 7 

/ Yl '^i-M^iA^^ ^) - 0' 0))df-^-^dijY—dxdT 

J to Jxo i—Q P+ 

<C / 5^||9^(a,,(-,F,a)-a,,(-,0,0))|| 



(210) 



L2[xo,L] 



42 



Using Cauchy-Schwarz inequality, we have 



rto+T PL'-—* -dZ 

/ Ci-Maij{x, y, a) - aij{x, 0, 0))d'^-^-'dijY^dxdr 

J to J XO l—Q P+ 



<c 



to+T 

to 

1/2 



A;-3 



$,_4(F,r;F,a) + 5^ 

1/2 



1=0 



da 


2" 


dt^ 





1/2 



X 



(211) 



^l'^,{Y, r; Y, a)^',^_\{Y - Y, r; Y, a)dT 

rto+T 

J to 

1/2 



^1^^,{Y, r; Y, a)e,ZiY, r; Y, a)^l^^,{Y - Y, r; Y, a)dT 



+ C 



to+T 



■fe-3 



d^a 



dt^ 



<^'J^,{Y, r; Y, a)<i>l^^,{Y - Y, r; Y, a)dr 
F, r; F, r; F, aj^^/.'^lF - Y, r; F, ajdr 



.1/2 



'to 

J to 

-to+T 

.V2 



+ C / $;/_2(F,r;F,a)$fc_4(F-F,T;F,a)(iT 

•/to 

J 1=0 



d^a 



dt^ 



^]^^^{Y,T-Y,a)dT. 



Similarly, we can estimate all other terms. In summary, we have the estimate as follows: 



$fc_4(F - F, to + T; F, a) + 4-4(F - F, to + T; F, a) - 4-4(F - F, to; F, a) 



<|.fe_4(l^ - *o; Y, a) + C|| II (F, a) II II 



d^a 



dt^ 



dr 



(212) 



+ C\ 

+ C\ 
+ C\ 

+ c\ 
+ c\ 



\{y,<^)\ 

\{Y,a)\ 



*o 1=0 

to+T 

^/^,{Y, r; F, a)|.;/_'4(F - F, r; F, ajdr 

to 
to+T 

dfe_4(F-F,r;F,a)(iT 

to 

(4-4(V' - F, to + T; F, a) - 4-4(V' - F, to; F, a)) 
(4-4(1^, to + T; F, a) - 4-4(1^, ^o; 1^, f^))- 
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Using the equation fl207p . we have 

A:-3 I 

E E - Y) Hi. <C^u-.{Y - F, r; a)(l + E E H^^^^^^H 

(213) 

+c||ii(v',^)ii|IE 



fc-3 I 



1=0 i=0 



2 > 
L2y 



i=0 



By the definitions of $o and ^o, there exists a C(ao) > such that 
33 + 31ao 



(214) 
and 

(215) 



64 



-^k-4{Y,to + T;Y,a) 



<^k-4iY, to + T- y, a) + C(«o)$fc-4(V' - F, to + T; F, a) 



128 



(4_4(F, to + T; r, a) - #fc_4(F to; F a)) 



<-(#fc_4(F, to + T; F, a) - #fc_4(F, to; Y, a)) 
o4 



+ C(^fc_4(r - F, to + T; F a) - ^fe_4(F - F, to; F ^)) 
where ao ^ (0, 1) is the constant in (11221) . Therefore, 
33 + 31^0 



(216) 



^, ^k-AiY, to + T; F cj) + — (^fc-4(F to + T; F ^) - ^fc-4(F to; F ^)) 
64 128 

<$fc_4(F, to + T; F t^) + C{ao)^u-^{Y - F, to + T; F o^) 

+ 7^(^fc-4(F, to + T; F a) - #fc_4(F, to; F ^)) 
64 



+ C(4-4(>^ - F, to + T; F fx) - ^k-A{Y - F, to; F ^)) 
Note that the Theorem [6] yields 



(217) 



<lfc_4(F, to + T; F ^x) + 7^(^fc-4(F, to + T; F ^^) - ^fc-4(F, to; F t^)) 

64 

<^^$.-4(r,to;Fa). 



Combining (12161) and (121 7p together yields 

^^±|i^<l,_4(F to + T; F ^) + Y^(^^-4(F to + T; F ^) 
(218) <li^|,^_^(y^ ^) + C{ao)^,^,{Y - F, to + T; F o) 
+ C(4-4(>^ - F, to + T; F cr) - 4-4(>^ - F, to; F t^)) 



_4(Fto;F^)) 
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Using fl200D . fl205D . and fl206D . we have 

33 + 31ao |,^_^^y^ ^ ^ _L(4_4(y, + T; y, a) - #fc_4(y, to; a)) 

^ l + 31ao ^^_^^y^ to; + C^k-4{Y - F, to; Y, a) 
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(219) 



+ C\ 
+ C\ 

+ c\ 
+ c\ 



\iY,CT) 

\iY,a) 
\iY,cr) 



a 



dr 



to+T k-3 
*o 1=0 



to 



<l>,_4(r-y,r;y,a)rfr 

(^fc_4(r, to + T; r, a) - 4-4(F to; Y, a)) 
{%-,{Y - F, to + T; r, a) - 4-4(>^ - F, to; F ^)) 



< 



1±^^^^_,{Y, to; F ^) + Ce<l>fc_4(F to; F 
+ Cel'fc_4(F,to;F^^) 
+ Ct ^k-4{Y,r-Y,a)dr 

J to 

+ Ce(4-4(F to + T; r, a) - 4-4(F to; F fx)). 
In the same way as (11971) . we have the following estimate 

l>fc_4(F t; F + (1 - Ce)(^fc_4(F t; F f^) - ^fe-4(F to; F ^)) 
<2<lfc_4(Fto;Ff^)- 

Substituting fl220|) into fl2T9|) and using fl205|l . and fl206D . one has 



(220) 



(221) ^ii|2^*.-.(y,to + r;y,.) < i±^*,_,(y,i„;y», 

o4 32 
if e is sufficiently small. As same as the proof for Corollary [9], it follows from (12211) that 
one has 



(222) 

where 

(223) 



l',_4(F, t; F s) + a\t) < C(<lfc_4(F 0; F (t) + (T2(0))e 



-2Af 



lni±^ ^ 17+15ao 
and a — 



2T 



2 + 30ao 



45 



Thus 



(224) 



fc-6 



5^||r(^,■)IU^[.,.i.]<C6V^^ 



1=0 



if we have ^k-^iY, 0; Y, a) < e^. This yields that 



(225) 



fc-6 

E 

i=0 



< Ce'e 



2-Xt 



Combining ( 1204 p . fl224p . and fl225p . one has fll63p . This finishes the proof of the Proposi- 
tion [101 □ 



Once one has the Proposition [TOl Theorem |2] follows from the standard continuation 
argument and local existence result. 



5. Linear Instability 

We consider the linear instability for the steady transonic shock solutions when E{xq) < 
—C for some positive constant C, where Xq G [0, L) is the shock location for the steady 
transonic shock solution. We consider the linearized problem (11210 . which can be written 
in the following form: 



f duY - ((c2(p+) - ul) d,Y) + {2u+dtY) + p+Y + E+d,Y = 0, 
for t > 0, < X < L, 



(226) 



d,Y 



_y2 

— -dxY —Y, at X = Xq, 

2^ I 2ua_ 



dxY = 0, dX X = L, 



where c^(p+) = p'{p+), u+ = Suppose that 



(227) 



E+{xo) < -C 
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for some positive constant C. We look for solutions for the problem fl226p of the form 
Y = e^^Z. Then 

r - ul)dlZ + (9.(c2(p+) - ul) - 2u^X - E^)Zx 



(228) 



- {X^ + 2Xd.xU+ + p)Z = 0, 



d.,Z 



C2(p+) - 

drZ = 0, at X = L. 



2m, 



for Xo < x < L, 



We use the shooting method. Given a parameter Z{xo) = a > 0, we consider the 
following initial value problem: 

f (p'(p+) - ul)dlZ + (9,.(p'(p+) - 4) - 2m+A - 

(A^ + 2XdxU+ + p)Z = 0, for x > Xq, 

2u, 



(229) 



Z{xo) = a>0,d,Z 



-2-(^o) ( + A 1 atx 



Xq. 



P'(P+) V 

Since > and c^(p+) > «+. Therefore, in view of f l227p . if A = 0, then Z^^xq) < 0. 
Therefore, there exists Li > xq such that Zx{x) < for xq < x < Li. 

If A = — ^(xo), then Zx{xq) > 0. then there exists L2 > xq such that Zx{x) > for 
Xq < X < L2. Let 

L = min{Li, L2}. 

By the continuous dependence of ODE with respect to the initial data and the param- 
eters, there exists a < A < — 22^(xo) such that the problem fl229p admits a solution 



Z satisfying Zx{L) = which is a solution of 



on [xo,i^]. This implies that the 



linearized problem can have exponentially growing solutions. 
This finishes the proof of Theorem [3l 



6. Appendix: Rauch- Taylor argument 

In this appendix, we give a detailed calculation for the estimate fll32p . which is mo- 
tivated by the estimate in [34j. For simplicity, we introduce the following characteristic 
variables for the equation in fll2ip : 

1 . 1 r 1 



(230) 



1 '•^ 

e = t-- 



c4 



-dx 



2 Jxo u+ + c(p+) 2 y^.„ u+ - c(p+) 



-dx. 



1 



-dx 



M+ + c(p+) 2 J^^ u+ - c(p+) 



-dx. 
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Obviously, C = if x = Xq. We also denote ( = (l if x = L. Then the equation in fll2ip 
will be 

(231) deeY - d(;(;Y + Md^Y + NY = 0, 

where 

Let Y = e^^Z, then one has 

(233) %Z - deeZ + {2k - M)a^Z + (fc^ - kM - N)Z = 0. 
If we Choose /c sufficiently large, then 

(234) 2fc - M > and - A;M - iV > for C e [0, Cl]- 

Multiplying the equation fl233p with di^Z and integrating on the domain fl = {{(,0)\9o + 
C < ^ < ^1 - C, < C < Co} with Co e [0, Cl], then we have 

(dcZfnc - 2dcZdeZng + (deZfric , f - kM - N)Z^ , 

— ^ ^ ^ -as + / 71(^0,3 

an 2 Jq^ 2 



{2k - M){d^Zy + ^ '-Z'dOdC = 0, 



(235) 

where {rig, ra^) is the unit outer normal on the boundary dQ. Therefore, we have 

(236) + t t~\2k - M){d,Z)\BX) + ^-^^^^^Z\QX)dBdC, 

JO JOo+C ^ 

^ r (M)!±(M)!„,o) + (il^M^i^(9.oM.. 

Jeo 2 2 

here we require Oi — 9q > 2(l- Integrating this inequality with respect to Co on [0, Cl] 
then we have 

(237) + f^'dCo r r \2k- M){d^zr{ex) + ^^^'^^^^z\ex)dedc 

Jo Jo Jdo+C ^ 
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If k is sufficiently large, then 



Jo Jeo+c. ^ 
Thus the estimate f l237p yields 



< I dc I ^ '—{ex)de. 



Jeo ^ ^ 

Transforming back to x and t coordinates, one has 

(238) ^dx f\{dtYf + {d^Yf + Y\t,x)dt<C f {{dtYf + {d^Yf + Y\t,xo)dt 

J xq J to Jo 

for some < to < < Note that 

(239) My,t)<C f\{dtYf + {d^Yf + Y\t,x)dx, 

J xo 

one has 

(240) r [{dtYf + {d^Yf]{t,xo)dt>C ['''\o{Y,t)dt-C Y'{t,xo)dt 
Jo J^-s Jo 

for T > suitably large and 6 G (0,T/4). This is exactly the estimate (I132p . 
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